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ABSTRACT 

Motivated by theoretical predictions that first stars were predominantly very massive, 
we investigate the physics of the transition from an early epoch dominated by massive 
Pop III stars to a later epoch dominated by familiar low-mass Pop II/I stars by 
means of a numerically-generated catalogue of dark matter halos coupled with a self- 
consistent treatment of chemical and radiative feedback. Depending on the strength 
of the chemical feedback, Pop III stars can contribute a substantial fraction (several 
percent) of the cosmic star formation activity even at moderate redshifts, z « 5. We 
find that the three z w 10 sources tentatively detected in NICMOS UDFs should be 
powered by Pop III stars, if these are massive; however, this scenario fails to reproduce 
the derived WMAP electron scattering optical depth. Instead, both the UDFs and 
WMAP constraints can be fulfilled if stars at any time form with a more standard, 
shghtly top-heavy, Larson IMF in the range IM© < Af* < lOOM©. 

Key words: galaxies: formation - cosmology: theory - cosmology: observations - 
intergalactic medium 



1 INTRODUCTION 

In the last few years, the universality of the stellar initial 
mass function (IMF) has been questioned by the results of 
■ theoretical studies, which consistently predict that the first 
] stars (hereafter Pop III stars) had characteristic masses of 
' 100-600 Mo (Omukai & Nishi 1998; Abel, Bryan & Norman 
2002; Bromm, Coppi & Larson 2002; Omukai & Palla 2003), 
about 100 times more massive than those observed today 
(Larson 2003). 

Indeed, the physical conditions in primordial star- 
forming regions appear to systematically favor the forma- 
tion of very massive stars. In particular, (i) the fragmenta- 
tion scale of metal-free clouds is typically 10"^ Mq (Abel, 
Bryan & Norman 2002; Bromm, Coppi & Larson 2002) (ii) 
because of the absence of dust grains the radiative feedback 
from the forming star is not strong enough to halt further 
gas accretion (Omukai & Palla 2003). (iii) since the accre- 
tion rate is as large as 10~''-10~^ Afoyr"^, the star grows up 
to ^ 100 Mq within its lifetime (Stabler, Palla & Salpeter 
1986; Omukai & Nishi 1998; Ripamonti et al. 2002). 

In spite of the significant progresses made by numerical 
simulations and semi-analytic models, many aspects of pri- 
mordial star formation remain to be fully understood, the 
more important of which are the radiative feedback effects 
which operate during protostellar accretion and which are 



likely to ultimately set the final stellar mass. By modeling 
the structure of the accretion flow and the evolution of the 
protostar. Tan & McKee (2004) have recently shown that 
radiative feedback becomes dynamically significant at pro- 
tostellar masses « 30 Mq , and are likely to constrain the 
mass of first stars in the range 100 — 300 Mq. 

On the other hand, observations of present-day stellar 
populations show that stars form according to a Salpeter 
IMF with a characteristic mass of 1 Mq, below which the 
IMF flattens. Thus, unless the current picture of primordial 
star formation is lacking in some fundamental physical pro- 
cess, a transition in the modes of star formation must have 
occurred during cosmic evolution. 

The physical processes responsible for this transition, as 
well as its cosmological consequences, are currently the sub- 
ject of systematic investigation. In particular, following the 
early study by Yoshii & Sabano (1980), it has been shown 
that a key element driving this transition is the metallicity 
of the star-forming gas (Omukai 2000; Bromm et al. 2001; 
Schneider et al. 2002, 2003). The fragmentation properties of 
the collapsing clouds change as the mean metallicity of the 
gas increases above a critical threshold, Za = 10~^^^ Zq 
(Schneider et al. 2002, 2003). The fragmentation of clouds 
with Z < Zci proceeds only down to relatively massive 
(> 100 Mq) cores, whereas in clouds with Z > Za lower 
mass fragments can be formed. Within the critical metallic- 
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ity range, low-mass cores can form if a sufficient amount of 
metals are depleted onto solid dust grains, which provide an 
efficient cooling channel in the high density regime (Schnei- 
der et al. 2003; Omukai et al. 2005). 

On cosmic scales, the transition between an epoch domi- 
nated by the formation of massive Pop III stars to an epoch 
dominated by the formation of ordinary Pop II/I stars is 
controlled by the strength of the so-called chemical feed- 
back: the explosions of the first massive supernovae starts to 
mctal-onrich the gas out of which subsequent stellar gener- 
ations form. In particular, in regions of the Universe where 
the metallicity exceeds Zct, low-mass stars can form. Chem- 
ical feedback sclf-propagatcs through the metal enrichment 
process and therefore depends on many poorly constrained 
parameters, such as the stax formation efficiency at high red- 
shifts, the number of supernovae per stellar mass formed, 
the efficiency of metal ejection and mixing. In particular, 
the shape of the primordial IMF governs the number of Pop 
III stars with masses in the pair-instability supernova range 
(PISN, Hcgcr & Woosley 2002) 140 Mq - 260 Mq (which 
contribute to the metal enrichment) relative to that of Pop 
III stars which end up as massive black holes (which instead 
do not contribute to the metal enrichment as all the heavy 
elements are captured by the black hole). 

Using a semi-analytic approach, Scannapieco, Schnei- 
der & Ferrara (2003) have modeled the probability that a 
collapsing halo has been impacted by metal-enriched out- 
flows from neighboring halos. Exploring different chemical 
feedback strengths, they concluded that the transition be- 
tween the Pop III and Pop II/I dominated star formation 
epochs is essentially controlled by the spatial distribution 
of metals. Because metal enrichment is generally inefficient 
(the filling factor of metals is always found to be < 1), the 
transition epoch is extended in time, with coeval Pop III 
and Pop II/I star formation episodes, and Pop III stars still 
contributing to the global star formation rate at redshifts 
z ^ 5. These findings have very important implications for 
the detection of Pop III stars through their strong Lya emis- 
sion (Scannapieco et al. 2003) as well as the corresponding 
pair- instability supernova rate (Scannapieco et al. 2005). 

In this work, wo investigate the cosmic transition be- 
tween Pop III and Pop II/I stars following a complemen- 
tary approach. We use the code PINOCCHIO* developed 
by Monaco, Theuns & Taffoni (2002) and Taffoni, Monaco & 
Theuns (2002) to produce random catalogues of dark matter 
halos at different redshifts. We then propagate metal enrich- 
ment along the hierarchy of mergers of star-forming galaxies 
to predict the number and mass distribution of metal-free 
objects, which can host Pop III stars, at each redshift. We 
then couple these findings to the star formation history (of 
Pop III and Pop II/I stars) and to the history of cosmic 
reionization and quantify their interplay through radiative 
feedback effects which can contribute to the extinction of 
Pop III stars. 

A similar approach has been recently undertaken by 
Purlanetto & Loeb (2005) with the aim of assessing the plau- 
sibility of the double reionization scenario suggested by Cen 
(2003) as a way to meet both the constraints on the gas 
neutral fraction implied by observations of Gunn-Peterson 

* http://www.daut.univ.trieste.it/pinocchio. 



troughs in the spectra of distant quasars (Fan et al. 2002, 
2003) and on the gas ionized fraction derived by the high op- 
tical depth to Thomson scattering measured by the WMAP 
satellite (Spergel et al. 2003; Kogut et al. 2003). 

Here we instead focus on the extinction epoch of Pop 
III stars which we believe to be regulated by feedback both 
of radiative and chemical types. We investigate the impli- 
cations that complementary observations on galaxy number 
counts at redshifts 2 « 10 (Bouwens et al. 2005) and on 
the history of cosmic reionization (Kogut et al. 2003) set 
on the transition epoch between Pop III and Pop II/I stars 
and on the nature of Pop III stars as sources of UV photons 
contributing to cosmic reionization. 

In particular, the detection of galaxies at the highest 
redshifts is challenging because of the redshifting of UV 
light into the infrared and the well-known limitations of 
current IR instruments. In spite of these difficulties, the 
search for high redshift galaxies, based on the dropout tech- 
nique, has provided the first tentative detection of galaxies 
at redshifts beyond z « 7 (Bouwens et al. 2004), stimu- 
lating a re-analysis of the small amount of deep IR data. 
Bouwens et al. (2005) looked at the prevalence of galax;ies at 
z ~ 8 — 12 by applying the dropout technique to the wide va- 
riety of deep fields that have been imaged with the Near In- 
frared Camera and Multi-Object Spectrometer (NICMOS) 
on board of the Hubble Space Telescope. Using an appro- 
priate selection criterion for high-z sources, Bouwens et al. 
(2005) concluded that the number of z ~ 10 sources must 
be three or fewer. When combined with the high optical 
depth to Thomson scattering inferred from the Wilkinson 
Microwave Anisotropy Probe (WMAP) satellite (Kogut et al 
2003), these data already provide a demanding benchmark 
for theoretical models. 

The paper is organized as follows. In Section 2 we de- 
scribe the different feedback effects that control the transi- 
tion from a Pop III to a Pop II dominated star formation 
epoch. We illustrate the adopted chemical feedback model, 
predicting the rmmbcr and mass function of halos which 
can host Pop III stars at different redshifts. In Section 3, we 
couple these predictions to a star formation and reioniza- 
tion model. In Section 4 we summarize the observations of 
galaxy counts at redshifts z = 10 and discuss how these can 
be used to set important constraints on the model. Finally, 
in Section 5 we present the main results and in Section 6 we 
discuss their implications. 

Throughout the paper, we adopt a ACDM cosmological 
model with parameters = 0.3, f^A ~ 0.7, h — 0.7, fis = 
0.04, n = 1 and as = 0.9, which lie within the experimental 
errorbars of WMAP experiment (Spergel et al. 2003). We 

also use the AB magnitude system^. 



2 FEEDBACK-REGULATED TRANSITION 

At least three types of feedback effects can be relevant to the 
extinction of the first massive stars: (i) radiative feedback 
caused by H2 photo-dissociating radiation emitted by the 

t AB magnitudes are defined as AB= —2.5 log(Fi,Q ) — 48.6, where 
is the spectral energy density within a given passband in units 
of erg cm~^ Hz~^ 
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Figure 1. Properties of the reference PINOCCHIO run. The hor- 
izontal line shows the mass of the smallest resolved halo and the 
line with points indicates the mass of the largest halo in the box 
at each simulation output. The solid line shows the evolution of 
the minimum halo mass allowed to form stars, corresponding to 
a virial temperature T^i^ = 10^ K. The short- and long-dashed 
lines represent the evolution of the filtering mass for two assumed 
reionization histories where reionization is complete at redshifts 
10 and f» 16, respectively (see Section 2.2). 



first stellar sources, which can iniiibit gas cooling and star 
formation inside small halos (ii) chemical feedback due to 
the release of metals and dust in the first Pop III supernova 
explosions, which pollute the gas out of which subsequent 
generations of Pop II stars form (iii) radiative feedback due 
to photo- heating filtering caused by cosmic reionization, i.e. 
the increased temperature of the cosmic gas suppresses the 
formation of galaxies below the Jeans mass. A thorough dis- 
cussion of each of these processes is beyond the scope of this 
paper and can be found in the recent extensive review by 
Ciardi & Ferrara (2005). 

As a first step, we build up the hierarchy of merg- 
ers of dark halos hosting the star formation sites. We use 
PINOCCHIO code (Monaco, Theuns & Taffoni 2002; Taf- 
foni, Monaco & Theuns 2002) to generate random catalogues 
of dark matter halos in a hierarchical universe. 

Our reference run has 256^ particles within a volume of 
comoving side 6 Mpc so that the smallest resolved dark 
matter halos have total mass 1.5 x lO'^ Mq (in the merger 
trees only halos with at least 10 particles are considered). 
Fig. 1 summarizes the details of our reference PINOCCHIO 
run and shows the redshift evolution of the largest and small- 
est mass present in the simulation box and of the minimum 
halo mass that is allowed to form stars. The latter corre- 
sponds, at each redshift, to halos with virial temperatures 



= lO'^M, 



-3/2 



(1) 



SO that stars do not form in halos which rely only on H2 
cooling (radiative feedback). The chosen resolution ensures 



Figure 2. Comparison between the predicted mass function of 
halos within the simulation box (points with errorbars) and the 
analytic predictions of Press & Schechter (1974, solid lines) and 
Sheth & Tormen (1999 dashed) at redshifts 30, 20, 10 and 5 (from 
left to right). Data at redshift 5 have been shifted upward by 
a factor 0.3 for clarity. The errorbars are Poisson errors on the 
counts in each mass bin. 



that we can follow the formation of stars in the lowest mass 
halos at the highest redshift. Due to radiative feedback, the 
first star-forming halos have a total mass « 2 x lO^M© and 
start to form at z ~ 30. This resolution is appropriate to 
follow the transition from a Pop III to a Pop II star forma- 
tion mode around redshift z « 10, which is the main scope 
of this analysis. 

In Fig. 2 we compare the resulting mass function of ex- 
isting halos at redshifts 30, 20, 10 and 5 in the simulation 
box, to the analytic predictions of Press & Schechter (1974) 
and Sheth & Tormen (1999). The errorbars are Poisson er- 
rors on the counts in each mass-bin. The agreement is very 
good down to redshift 0^5 showing that the chosen reso- 
lution ensures completeness in the simulation volume in the 
redshift range relevant to our analysis. 



2.1 Chemical feedback 

At each redshift, existing halos which are allowed to form 
stars (Tvir ^ lO** K) are classified as hosting Pop II (Pop 
III) stars depending on whether the halo itself or any of its 
progenitors have (have not) already experienced an episode 
of star formation, leading to the formation of SN progenitors 
which then metal-enrich the gas. Note that even a single 200 
A'Iq PISN can enrich the gas within a halo with mass Mmin 
to a metallicity level well above the critical value as. 



0.45M. 



{nB/nM)M„ 



0.02Z, 



(2) 



where, following Heger & Woosley (2002), we have as- 
sumed that the mass of metals released in a PISN explosion 
is 45% of the initial progenitor mass. 
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This way of taking into account chemical feedback is 
very conservative as (i) we are assuming that all metals syn- 
thesized during a star formation episode remain confined 
(or are re-accreted) within the host halo and/or (ii) even if 
a fraction of metals is ejected out of the host halo, the metals 
that are kept within the potential well are sufficient to pol- 
lute the gas within the host system or within any unpolluted 
halo that the system will later merge with, to a metallicity 
level larger than Zcv On the other hand, we are neglecting 
metal pollution from neighboring galaxies through enriched- 
winds (Scannapieco et al. 2003; Sigward, Ferrara & Scanna- 
pieco 2005). 

Fig. 3 shows the fraction of halos which can host Pop III 
stars (hereafter Pop III halos) at different redshifts as a func- 
tion of their mass. The three different curves correspond to 
different assumptions on the shape of the primordial IMF. In 
particular, we randomly classify Pop III halos as SN-forming 
or BH-forming depending on whether they host stars with 
masses in the PISN range or outside this range. SN-forming 
halos are metal-enriched and pollute their descendants. BH- 
forming halos are not enriched in metals and therefore do 
not propagate chemical feedback along their subsequent hi- 
erarchy of mergers. The fraction of SN-forming halos, /sn, 
is therefore a free parameter that controls the strength of 
chemical feedback. 

The three curves shown in each panel of Fig. 3 corre- 
spond (from bottom to top) to /sn = 1,0.5,0.1. Each curve 
represents an average over 50 realizations of the random 
classification procedure and the errors are Poisson errors on 
the counts in each mass bin. 

When /an = 1, all Pop III halos are SN-forming and 
chemical feedback is strong. The total fraction of Pop III 
halos decreases with time, being « 0.25 at 2; = 20 and ~ 0.1 
at 2 = 5; at all redshifts, Pop III stars are confined to form in 
the smallest halo mass bins (close to the minimum mass to 
form stars shown as the vertical dotted line in each panel). 
For /sn = 0.5 (50% of Pop III halos host SN), the number 
of Pop III halos at each redshift is larger than the previous 
case (~ 0.5 at z = 20 and « 0.2 at z « 5) but Pop III ha- 
los are still confined to the smallest mass bins ( < 10^ M0). 
Finally, when /sn = 0.1 the strength of chemical feedback is 
significantly reduced (as only 10% of Pop III halos at each 
redshift host SN) and therefore the fraction of Pop III ha- 
los is dominant at all redshifts, being « 0.8 at 2; = 20 and 
« 0.65 at 2 = 5, not negligible even in halos with the largest 
mass present in the simulation box at each redshift. Here- 
after we only consider /sn = 0.1 and = 1 as representative 
of weak and strong chemical feedback models. 

Despite the strength of chemical feedback. Pop III halos 
continue to form stars even at redshifts as small as 2^5, 
in agreement with the results of Scannapieco et al. (2003). 
Thus, independently of whether it preferentially propagates 
through galaxy mergers or through metal-enriched winds, 
chemical feedback always leads to a smooth transition be- 
tween a Pop HI and a Pop II dominated star formation 
epoch. Note that the above results depend only on the pa- 
rameter /sn and on the hierarchy of galaxy mergers, and are 
totally independent of the star formation efficiency or of the 
value of critical metallicity for the transition in fragmenta- 
tion scales of collapsing pre-stellar clouds. 




lOa iQS io'» 10" 10= 10" lO'o 10" 
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Figure 3. Fraction of star forming halos which, at a given 
redshift, hosts Pop III stars as a function of their total mass. 
The different panels correspond to different output redshifts, 
z = 20, 15, 10 and 5. In each panel, the three curves correspond 
to /sn = 0.1,0.5, 1 (from top to bottom) and the vertical dotted 
line shows the mass that at each considered redshift corresponds 
to the minimum halo mass to form stars, M = Mmin (see text). 

2.2 Photoionization 

So far, we have taken into account radiative feedback only 
through the choice of the minimum halo mass that can 
host stars. Namely, we have assumed that mini-halos with 
virial temperatures less than 10* K do not significantly con- 
tribute to the cosmic star formation history. However, ra- 
diative feedback does affect larger halos, particularly when 
the universe is approaching the epoch of cosmic reioniza- 
tion. Indeed, the increase in temperature of cosmic gas dur- 
ing reionization will suppress the formation of small galaxies 
with mass below the Jeans mass. To quantify this effect, we 
follow the method proposed by Gnedin (2000) and adopt the 
filtering mass as the characteristic scale below which the gas 
fraction in star forming objects, hence the fraction of stel- 
lar mass, is reduced with respect to the universal value. In 
particular, before reionization, the average baryon fraction 
(the ratio of the mass in baryons over the total halo mass) 
has the universal value Q,b/^m while after reionization it 
can be approximated by the fitting formula (Gnedin 2000), 

r ^ 0.3/ , . 

[1 + (2V3 _ l)Mc/M]3 ' ^ ' 

where the characteristic mass. Ale, is the total mass of 
objects which, on average, retain 50% of their universal 
baryon mass, 0.5(^Ib /^m)Mc ■ This characteristic mass is 
very well approximated by the linear-theory filtering mass, 
Mc « Mf, with 

M^/«=3 rda'M'/'iT,a')[l^i^)'^% (4) 

"■Jo " 
where the Jeans mass is given by 

Mj{T) = 2.56 X 10^ MQ{nMh'^r^^^iT/fif^\l + 2)-^/^(5) 
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Table 1. Properties of the three models considered: fsn is the strength of chemical feedback (the fraction of Pop III halos which hosts 
SN and therefore contributes to metal enrichment), columns 3-5 (6-8) show the initial stellar metallicity, IMF type and Log of Hl-iouizing 
photons per baryons in stars for Pop III (Pop II) stars. By Top-Heavy IMF we mean that only stars with masses ^ 100 Mq are considered. 
Finally, the last two columns show the minimum and maximum values for (/esc, /*) allowed by the observed NICMOS-UDFs counts at 
redshift 2 = 10 (see text). 









Pop III 






Pop II 








Model 


fsn 


Z/Zq 


IMF 




Z/Zq 


IMF 




[feso f*]min 


[fesci f-klmax 


A 


0.1 





Top-Heavy 


4.862 


0.2 


Salpeter 


3.602 


[0,0.9 X 10-3] 


[0.5, 1.65 X 10-3] 


B 


1 





Top-Heavy 


4.862 


0.2 


Salpeter 


3.602 


[0,1.15 X 10-2] 


[0.5,2.15 X 10-2] 


C 
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with T and /j the temperature and molecular weight of the 
cosmic gas. 

Here we assume the temperature of ionized regions to 
be 10'' K and we modulate the increase in temperature of the 
cosmic gas due to reionization with the volume filling factor 
of ionized regions at the corresponding redshift Qhii(2;) (see 
Section 5.1). 

Thus, the characteristic mass depends on the thermal 
history of the gas and it accounts for the finite time required 
for pressure to infiuence the gas distribution. The redshift 
evolution of the filtering mass is shown in Fig. 1 for a model 
where reionization occurs at redshift Zrei ~ 10 (short-dashed 
line) and ~ 16 (long-dashed line). The filtering mass steeply 
increases with time and, as expected, its value in the red- 
shift range 5 ^ z ^20 depends strongly on the reionization 
history. At redshift z « 10, the filtering mass is « 10® M© 
for the late reionization case, and it is a factor ten larger for 
the early reionization case. 

As can be inferred from Fig. 1 and 3, following reion- 
ization, photoionization feedback can significantly suppress 
the formation of Pop III stars since halos with masses 
10* Mq ^ M ^ IO^Mq have evaporated at least 50% of 
their mass. 



3 REIONIZATION HISTORY 

Following Barkana & Loeb (2001), we compute the redshift 
evolution of the filling factor of HII regions as. 



Qhii(^;) 



«^ z 



dz 



1 dn 
rK~dt 



(6) 



where ri^ = XHn% and Ug are the present-day number den- 
sities of hydrogen and baryons {Xh ~ 0.76 is the hydrogen 
mass fraction), and duj/dt is the production rate of ionizing 
photons that we can express as, 



drij _ 
dt 



(7) 



where fesc is the escape fraction, i. e. the fraction of ionizing 
photons that is able to escape the host galaxy, Nj is the 
time-averaged number of ionizing photons emitted per unit 
stellar mass formed, and p« is the comoving star formation 
density, 



dMM-^{M,z) 
dM^ ' 



(8) 



where /« is the star formation efficiency, /a is the baryon 
fraction given by eq. (3), M-min is the minimum halo mass 
to form stars (see eq. 1) and n(M, z) is the number density 
of halos (Pop II and Pop III) at redshift z. The parameter 
fesc also affects the relative strength of the stellar continuum 
and the Lya line and nebular luminosity of each galaxy. The 
latter contributions are particularly important for massive 
Pop III stars. 

Finally, the function F{z' , z) takes into account the 
effect of recombinations. Assuming a time-independent 
volume-averaged clumping factor C, common to all HII re- 
gions, we can write, 



F{z',z) = 



C[fiz') - f{z)], 



and 



f{z) = J{l + z)^ + 



(9) 



(10) 



is the hydrogen recombina- 



where as = 2.6 x 10~ cm 

tion rate. 

For a given reionization history, it is straightforward to 
compute the optical depth to electron scattering as, 



Te{z) 



r dz' 


dt 






Jo 


"d^ 



aTne(z'), 



(11) 



where ctt = 6.65 x 10~^^cm^ is the Thomson cross Sec- 
tion and ne{z) is the number density of free electrons at 
redshift z. In this paper we only consider the evolution of 
the filling factor of HII regions (Qhii). However, for the 
purpose of computing the optical depth to electron scat- 
tering, we assume that Qneii = Qhii and we neglect the 
contribution of electrons coming from Helll regions (see also 
Haiman & Holder 2003). Under these assumptions we can 
write ne{z) = Qmi{z)n%{l + z)^ where Qhu{z) is given by 
eq.(6). 



4 GALAXY COUNTS AT REDSHIFT 10 

Bouwens et al. (2005) looked at the prevalence of galaxies 
at 8 — 12 by applying the dropout technique to the wide 
variety of deep FllOW- and F160W-band (hereafter Jno 
and Hum, respectively) fields that have been imaged with 
the Near Infrared Camera and Multi-Object Spectrometer 
(NICMOS). The principal data set is constituted by two 1.3 
arcmin^ deep NICMOS parallels taken with the Advanced 
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Camera for Surveys (ACS) Hubble Ultra Deep Field (UDF; 
each has about 160 orbits data). Complementary shallower 
fields possessing similar Jno + Hiao imaging have been an- 
alyzed Eis well. The 5a limiting magnitude for the UDFs 
is ~ 28.6 in Jno and ~ 28.5 in iJieo (0.6" aperture). Since 
these limiting magnitudes are reached only in the deepest re- 
gions of the two NICMOS UDFs, we adopt the conservative 
limiting magnitude Hieo = 28 (Bouwens, private communi- 
cation) . 

The primary selection criterion for liigli-2: sources is 
Jiio — Hieo > 1.8. Using this criterion, Bouwens et al. (2005) 
found eleven sources. Eight of these are ruled out as cred- 
ible ~ 10 sources, either as a result of detection (> 2a) 
blue-ward of Jno or because of their colors red-ward of the 
break (i/ieo — K ~ 1.5). The nature of the three remain- 
ing sources could not be assessed from the data, but this 
number appears consistent with the expected contamination 
from low-redshift interlopers. Hence, Bouwens et al. (2005) 
concluded that the actual number of 2: ~ 10 sources in the 
NICMOS parallel fields must be three or fewer. 

Adopting the same selection criterion as in Bouwens et 
al. (2005), we compute the number of 2 « 10 galaxies de- 
tectable in the NICMOS UDFs predicted by our models. 
Assuming that the three candidates selected by Bouwens et 
al. (2005) are indeed z k, IQ galaxies, we can set interesting 
constraints on the models. In particular, for a given chem- 
ical feedback strength /sn, we compute the number density 
of Pop HI and Pop II halos at each redshift. We then as- 
sume a universal value for the star formation efficiency /* 
and escape fraction f^sc of ionizing photons and we charac- 
terize each stellar population (Pop III and Pop II/I stars) 
by assuming a specific template emission spectrum. 

In the following, we consider two models, A and B, 
whose details are given in Table 1. We assume two repre- 
sentative chemical feedback strengths (/sn =0.1 for model 
A and /sn = 1 for model B). In both cases, the emission 
properties of Pop II stars have been computed using a sim- 
ple stellar population model taken from the GALAXEV li- 
brary (Bruzual & Chariot 2003). Pop II stars have an initial 
metallicity oi Z = 0.2 Zq and are assumed to form according 
to a Salpeter IMF (with masses in the range 0.1 Mq - 100 
Mq). Pop III stars are assumed to be very massive (with 
masses larger than 100 Mq) and their emission properties 
are computed using the stellar spectra of Schaerer (2002), 
including the nebular continuum emission found to be very 
important for stars with strong ionizing fluxes'^. 

For each model, we then compute the number of sources 
per unit solid angle with observed flux in the range Fi,^ and 
Fv„ + dF^„ as 



dN 



dQdF, 



"0 



iK,,zo)= rdz{^)niK„z), (12) 



where dV/dQ,dz is the comoving volume element per unit 
solid angle and redshift, n(F^o,a) is the comoving num- 



t Since the emission spectrum of very massive Z=0 stars is found 
to be independent of the stellar mass for stars with masses 
100 Mq (Bromm, Kudritzki & Loeb 2001; Schaerer 2002), we are 
not forced to make any assumption about the shape of the IMF 
of Pop III stars as long as these are assumed to be very massive. 



ber of objects at redshift z with observed fiux in the range 
[Jlo, Fi,„ + dF^g], given by 

n(F.„ ,z) = j^ dz' -^(F.„ , .') ^(M, .'). (13) 

Here dPn/ dMdz is the halo formation rate, and the factor 
dM/dFva converts the number density per unit mass in unit 
fiux interval. 

The average flux Fvq from a halo of mass M at redshift 
z that has formed at 2' ^ is. 



A(nB/OAf)M 

AixiluodLiz)^ 



di^liu,t,,,,)e-^-ff^''°''°-^\{U) 



where Z(!^, f^_^/) is the template luminosity per unit solar 
mass for a population of age t^^^i (the time elapsed between 
the redshift z' and z), dhi^z) is the luminosity distance, Afo 
is the instrumental bandwidth, Vmin and Umax are the rest- 
frame frequencies corresponding to the observed ones. Fi- 
nally, Tef f is the intergalactic medium (IGM) effective opti- 
cal depth at fo between redshift zo and z (see Sect. 2.2 of 
Salvaterra & Ferrara 2003) 



5 RESULTS 

Wc compute the expected number of z ~ 10 galaxies de- 
tectable in the NICMOS UDFs using eqs. 12-14 and apply- 
ing the same selection criterion of Bouwens et al. (2005) 
(Jno — Hieo > 1-8, integrated up to the magnitude limit of 
the survey, Hieo = 28). Assuming that only the three most 
luminous sources should be above the detection threshold, 
we derive the maximum allowed value for the star forma- 
tion efficiency, /,, corresponding to a given escape fraction 
of ionizing photons. 

The results obtained for models A and B are shown in 
Fig. 4 in terms of comoving star formation density in Pop HI 
and Pop II stars as a function of redshift. Model C will be 
discussed in section 5.1. 

In each panel, vertically (horizontally) shaded regions 
represent the range of star formation histories for Pop III 
(Pop II) stars allowed by the NICMOS-UDFs counts when 
a value ^ /esc ^ 0.5 in assumed (see Table 1 for the corre- 
sponding values of /*). For both models and all the explored 
values of fesc, the three most luminous objects in NICMOS 
UDFs are always represented by galaxies hosting Pop III 
stars, as their strong Lya and nebular emission (Schaerer 
2002) is dominant. As a consequence, /« is restricted to 
relatively small values, which increase with increasing fesc- 
Model B, being characterized by the strongest chemical feed- 
back, and therefore by the smallest fraction of halos hosting 
Pop III stars at redshift ^ 10 (see also Fig. 3), has larger 
allowed star formation efliciencies than model A. 

Fig. 4 shows that in model A chemical feedback is very 
inefficient (/s„ — 0.1) and Pop HI stars dominate the cosmic 
star formation history down to redshift z ^ 5. Conversely, 
model B shows that when chemical feedback is maximally 
efficient, the star formation density in Pop III stars is sub- 
dominant at all but the highest redshifts, but Pop III stars 
still make a significant contribution to the cosmic star for- 
mation history down to redshift 5. 

Because of the moderate star formation efficiencies im- 
plied by the observed NICMOS UDFs counts, radiative feed- 
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Figure 4. Star formation densities as a function of redshift in 
model A (upper panel), model B (middle panel) and model C 
(lower panel, see Table 1). The vertically (horizontally) shaded 
regions represent the range of star formation densities in Pop III 
(Pop II) stars allowed by the observed NICMOS-UDFs counts 
at redshift z = 10 (see text). For all models, the IGM clumping 
factor is assumed to be C = 10. 



Figure 5. Isocontours corresponding to reionization redshifts 6, 
8, 10, and 12 for model A in the plane (C, fesc), where C is the 
clumping factor and fesc the escape fraction. The dashed region 
indicates the parameter space which leads to an optical depth to 
electron scattering within the range Te = 0.16 ± 0.04 observed by 
WMAP. The dashed line corresponds to the central value. 



back is negligible for the most massive galaxies. To evaluate 
its impact on the predicted star formation densities shown 
in Fig. 4, we have assumed a clumping factor C — 10. For 
^ fesc ^ 0.5, the corresponding reionization redshifts vary 
in the range Zrci ^ 10 for model A, and in the range 
^5 Zrci ^ 14 for model B. 

In model A, radiative feedback effect on the overall 
Pop III star formation density manifests itself as a progres- 
sive decrement of the upper limit (which corresponds to the 
largest fesc and /*) of the shaded region for z ^ 10. In model 
B, where the star formation efficiencies implied by the ob- 
served counts are larger and, at each redshift. Pop III stars 
are confined to form only in the smallest halos, radiative 
feedback is more significant; around redshift 10, the lines 
corresponding to the maximum and minimum Pop III star 
formation efficiencies cross each other (the highest the star 
formation efficiency and escape fraction, the strongest the 
radiative feedback). 

Overall, the observed galaxy counts at redshift 2 « 10 
favor a cosmic star formation history where the transition 
from Pop III to Pop II stars is governed by strong chemi- 
cal feedback, with all Pop III halos efficiently propagating 
metal enrichment to their descendants along the hierarchy 
of galaxy mergers. In fact, in less efficient chemical feedback 
scenarios {fs„ = 0.1), the predicted cosmic star formation 
density is dominated by Pop III stars down to redshift z ~ 5. 

5.1 Implications for reionization 

Stronger constraints on models A and B come from their 
corresponding reionization histories. Figs. 5 and 6 show 
isocontours of redshifts of reionization as a function of the 



clumping factor and the escape fraction (the star forma- 
tion efficiency is set by the observed counts at z « 10). 
We define the redshift of reionization, Zrd, as the redshift 
when the filling factor of ionized regions, given by eq. (6) is 
Qmi{zTci) ~ 1- The shaded regions indicate the parameter 
space which leads to an optical depth to electron scattering 
within the range = 0.16 ± 0.04 observed by the WMAP 
satellite (Spergel et al. 2003; Kogut et al. 2003), with the 
dashed line corresponding to the central value. 

As it can be inferred from the figures, both models 
lead to reionization histories which are in poor agreement 
with observations. Indeed, large escape fractions and mod- 
erate clumping factors are required to reproduce the ob- 
served range of optical depths. Independent observations of 
the escape fraction give, with few exceptions, values in the 
range fesc ^ 0.15 (Ciardi & Ferrara 2005 and references 
therein) whereas numerical simulations (Gnedin & Ostriker 
1997; Iliev, Scannapieco & Shapiro 2005) and semi-analytic 
models (Madau, Haardt & Rees 1999) predict clumping fac- 
tors in the range 1 < C ^ 30. 

Fig. 5 shows that in model A, even assuming an ho- 
mogeneous IGM (C — 1), the escape fraction has to be 
fesc > 0.15 to be marginally consistent with the observa- 
tions. In model B, the agreement improves but still requires 
extreme assumptions on either the escape fraction or the 
clumping factor: when fesc ^ 0.15, the predicted optical 
depth lies within -la error from the observed value only if 
C < 8. 

Our analysis suggests that it is very hard to recon- 
cile the observed number counts at redshift z m 10 with 
the reionization history inferred by WMAP observations if 
Pop III stars are assumed to be all very massive, with masses 
^ 100 Mq . This result is somehow counter-intuitive because 
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t. 



Figure 6. Same as Fig. 5 but for model B and reionization red- 
shifts 8, 10, 12, 14 and 16. 



metal-free very massive stars have always believed to play 
an important role in an early epoch of cosmic reionization 
(Cen 2003; Haiman & Holder 2003; Sokasian et al. 2004, 
2004; Venkatesan, Tumlinson & ShuU 2003; Wyithe & Loeb 
2003; but see also Ciardi, Ferrara & White 2003, and Ricotti 
& Ostriker 2004 for alternative interpretations). 

Still, due to their strong Lya and nebular luminosity, 
Pop III stars with masses 100 Mq always dominate the 
predicted galaxy number counts at z ~ 10, even assum- 
ing the strongest chemical feedback model. The paucity of 
candidate objects observed in NICMOS UDFs implies that 
star formation in halos with T^ir ^ 10* K must have oc- 
curred very inefficiently. Even if all three objects identified 
by Bouwens et al. (2005) are truly z ~ 10 sources, the con- 
tribution of Pop III stars to cosmic reionization is strongly 
constrained. 

The constraints are particularly severe when chemical 
feedback is less efficient (model A), because Pop III stars 
can form at redshift z ~ 10 in relatively large, and therefore 
luminous, galaxies. Even if chemical feedback is assumed to 
be maximally efficient and Pop III stars are confined to form 
in the smallest halos, their contribution to cosmic reioniza- 
tion is in agreement with observations only for moderate 
(C < 8) clumping factors. 

Unless an additional type of feedback comes into play, 
the transition between Pop III to Pop II halos, regulated by 
chemical and radiative feedback, is very smooth. Thus, to 
relax the severe upper limits placed on the star formation 
efficiencies by NICMOS UDFs counts, the emission proper- 
ties of Pop III stars have to be modified. 

In particular, if Pop III stars are assumed to form 
according to a Salpeter IMF with masses in the range 
1 Mq ^ Mi, ^ 100 Mq, their Lya and nebular emission 
is strongly suppressed. In this limit, (see model C in Ta- 
ble 1), all Pop III halos contribute to chemical feedback, 
fsn = 1, because each burst of star formation leads to the 
formation of SN progenitors and NICMOS UDFs counts at 




t. 



Figure 7. Same as Figs. 4 and 5 but for model C and reionization 
rcdshifts 6, 8, 10, 12, 14 and 16. 

z « 10 are dominated by Pop II galaxies. As a consequence, 
the upper limit on the star formation efficiency set by the 
observed counts is /* = 0.1, much higher than those found 
for models A and B and independent of the escape fraction. 

The predicted evolution of the star formation density 
in Pop II and Pop III stars is shown in the lowest panel of 
Fig. 4. Similar to model B, Pop II star dominate the star 
formation history at all but the highest redshifts. Indeed, 
in our approach the transition between Pop III and Pop II 
dominated star formation epochs is controlled only by the 
strength of chemical feedback (/sn) and it is independent 
of the assumed star formation efficiencies (see section 2.1). 
However, contrary to model B, the star formation density 
depends on the escape fraction only indirectly, through the 
strength of radiative feedback. For this reason, the shaded 
regions representing the range of allowed star formation den- 
sities are not degenerate only after reionization which, for 
C = 10 and /esc = 0.5 occurs at redshift Zrd ~ 15 (for Pop 
II stars, which are preferentially hosted in large halos, this 
effect is negligible). 

The implications for the history of reionization are sum- 
marized in Fig. 7, which shows that model C is in much bet- 
ter agreement with observations than model A and B. The 
optical depth observed by WMAP can be reproduced with 
fesc ^ 0.15 even if the clumping factor is as high as C ^ 25. 

Thus, when the observed z ~ 10 counts are used to 
set an upper limit to the allowed star formation efficiency, 
Pop HI stars forming with a conventional IMF are more effi- 
cient sources of cosmic reionization than very massive stars. 
Note that this conclusion is independent on issues related 
to chemical and radiative feedback. Indeed, both models B 
and C assume maximally efficient chemical feedback and ra- 
diative feedback plays a minor role for model C because the 
predicted z « 10 counts are dominated by massive Pop II 
galaxies, which are largely unaffected by photoevaporation. 

Fig. 7 shows that the best-fit value of optical depth im- 
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Figure 8. Optical depths as a function of the characteristic stellar 
mass of Pop III stars, Mc, assumed to form with a Larson IMF 
(see text). The different shaded regions correspond to C = 1, 10, 
and 30. The emission properties of Pop III have been modified 
according to each assumed value of Mc and /esc = 0.1 has been 
adopted. For each C, the upper and lower limits to the shaded 
region corresponds, respectively, to 3 and 2: ~ 10 galaxies in 
NICMOS UDFs. 



plied by WMAP observations, Te = 0.16, can only be repro- 
duced by model C with /esc ^ 0.15 if C ^ 8. If larger values 
of C are to be favored or the observed optical depth will be 
corrected up-ward in future WMAP data releases, a revision 
of model C is to be considered. 

However, unless the NICMOS UDFs counts have under- 
estimated by a large factor the number of z « 10 sources, 
our results suggest that the only plausible revision of model 
C is to modify the primordial IMF by increasing the char- 
acteristic mass but allowing Pop III stars to form only in 
the conventional mass range, 1 Mq ^ K'h ^ 100 Mq. In 
Fig. 8 we show the predicted optical depth as a function of 
the Pop III characteristic stellar mass, Mc. We assume that 
the characteristic stellar mass enters in the definition of the 
stellar IMF as suggested by Larson (1998), 



dN 



dlogMv, 



l^+Mc) 



(15) 



which, for — leads to the standard Salpeter IMF. 

For a given Mc, we compute Te considering the proper 
emission rate of ionizing photons (taken from Table 1 of 
Ciardi & Ferrara 2005), assuming /esc = 0.1, three different 
clumping factors C = 1, 10, and 30, and that either all 3 or 
none of the candidates identified by Bouwens et al. (2005) in 
NICMOS UDFs represent z ~ 10 galaxies (upper and lower 
limits to each shaded region). Fig. 8 shows that the predicted 
optical depth depends more on the clumping factor than on 
the assumed value of Mc- If C = 10, the predicted Te is in 
agreement with the observed values even assuming Mc = 
0, consistent with the results shown in Fig. 7. If C = 30, 
WMAP and NICMOS UDFs data can be reconciled only 
if all three candidates are confirmed to be 2: ~ 10 galaxies 



and if the characteristic Pop III stellar mass is ^ 20 Mq. 
Finally, if the observed Tc will be corrected up-ward in future 
analyses of WMAP data, none of the models considered, 
independently of the assumed Mc, will be able to reconcile 
the observed reionization history with the counts ai z m 10, 
unless C < 10. 



6 DISCUSSION 

The tentative detections of 2 ~ 10 galaxies in NICMOS 
UDFs by Bouwens et al. (2005) provides tight constraints 
on the relative contributions of Pop III and Pop II stars to 
the galaxy luminosity function. When coupled to the reion- 
ization history of the universe implied by the WMAP satel- 
lite, this collection of data sets represents already a very 
demanding benchmark for galaxy evolution models. 

Stimulated by these data, we have followed the assem- 
bling history of high-redshift galaxies hosting the first stars 
in a hierarchical model for structure formation. Adopting 
simple prescriptions for chemical and radiative feedbacks 
along with available data, we have constrained the evolu- 
tion of the IMF of the first stars. In particular we find that: 

(i) Two types of feedbacks need to be considered: (ia) 
radiative feedback, which suppresses star formation in H2- 
cooling halos and the formation of low-mass galaxies be- 
cause of photoionization heating following reionization, and 
(ib) chemical feedback due to the progressive enrichment of 
(Pop II) star forming gas by the metals released by the first 
(Pop III) stellar explosions. 

(ii) If chemical feedback is weak {fsn ~ 0.1) Pop III stars 
can form in relatively large halos and dominate the evolution 
of the cosmic star formation density down to redshift z ^ 5. 
More realistic scenarios are found assuming /s„ = 1, where 
Pop III stars are confined to form only in the smallest halos 
at each redshift. 

(iii) The combined effect of radiative and chemical feed- 
backs fails to suppress the formation of massive Pop III stars 
at z < 10. 

(iv) If Pop III stars are assumed to form with a Top- 
Heavy IMF, that is with masses M* ^ 100 Mq, Pop III 
galaxy are always found to dominate the predicted z ^ 10 
NICMOS UDFs counts. Because of their strong Lya and 
nebular luminosities, very tight upper limits on the star 
formation efficiency in T^ir 10"* K halos are placed by 
the paucity of observed counts. The corresponding reioniza- 
tion histories fail to reproduce the observed WMAP optical 
depths unless uncomfortably large values for the escape frac- 
tion of ionizing photons and/or too small values for the IGM 
clumping factor are assumed. 

(v) Reionization constraints from WMAP and observed 
z^ 10 counts in NICMOS UDFs can be reconciled if Pop III 
stars form according to a Larson IMF with stellar masses 
in the conventional range 1 AIq Af* ^ 100 Mq and a 
characteristic mass Mc ^ if the clumping factor is 10, 
or Mc ^ 20 Mq if the clumping factor is as large as 30. In 
the latter case, at least 1 candidate identified by Bouwens 
et al. (2005) has to be a galaxy at z ~ 10. 

Our analysis is based on the assumption that chemi- 
cal feedback propagates along the hierarchy of galaxy merg- 
ers from progenitors to descendants rather than through 
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mctal-cnrichcd galaxy outflows (Scannapieco ot al. 2003). By 
adopting this simplified scheme, we are implicitly assuming 
that (i) even if a fraction of metals is ejected out of the host 
halo, the metals that are kept within the potential well are 
sufficient to pollute the gas within the host system or within 
any unpolluted halo that the system will later merge with, to 
a metallicity level larger than Zct and that (ii) metal pollu- 
tion from neighboring galaxies is negligible. In spite of these 
simplifying assumptions, the resulting Pop III and Pop II 
star formation histories arc consistent with the predictions 
of Scannapieco et al. (2003). In particular, even assuming a 
maximally efficient chemical feedback, the transition is al- 
ways predicted to be very smooth, with Pop III stars which 
contirmc to form in the smallest halos at each rcdshift down 
to z ^ 5. Interestingly enough, our conclnsion is indepen- 
dent of the primordial IMF and/or star formation efficiency. 
We show that models where all Pop III halos are able to 
metal pollute their descendants (/s„ — 1) are favored by 
observations. The smoothness of the cosmic IMF transition 
is therefore strictly linked to the inhomogeneous character 
of the metal enrichment process, independently of how it is 
implemented in the models. 

Even if we have adopted a very approximate descrip- 
tion of the thermal evolution of the cosmic gas, we find that 
radiative feedback following reionization can play an impor- 
tant role in the extinction of Pop III stars. In particular, if 
reionization occurs at redshifts « 15 (a « 10) the baryon 
fraction in 10^ M© halos is reduced to 70% (99%) of its uni- 
versal value (ns/fiM) at redshift 10 and to 10% (30%) of 
this value at redshift 5. Thus, the baryon content in the 
smallest halos where Pop III stars can form at low redshift 
is severely limited by radiative feedback. This is strongly 
related to the epoch of cosmic reionization which, in turn, 
depends on the degree of inhomogencitics which characterize 
the cosmic gas. We have parametrized this quantity through 
a constant clumping factor C and we have explored different 
values of C ranging from 1 to 50. However, it is likely that 
the degree of dumpiness of the cosmic gas varies with time, 
rising from C ~ 1 at high rcdshift to a few tens at Ri 10 
as structures go increasingly non-linear. 

Recently, Iliev et al. (2005) have studied the effect of 
IGM clumping at redshifts z ^ 10, providing a fit to the 
evolution of the clumping factor (see their Fig. 7) where 
C w 2 at redshift 30 and raises to C w 8 at redshift 10. 
Adopting this evolution for the clumping factor and fesc = 
0.1 we find that models A and B are unable to reproduce the 
optical depths inferred from WMAP data whereas in model 
C we find that reionization is complete at redshift 13.5 and 
that the resulting optical depth is Te = 0.1734 in very good 
agreement with the observed value. 

Thus, our analysis suggests that Pop III stars form- 
ing with masses in the same dynamical range as Pop II/I 
stars, but possibly with higher characteristic stellar mass, 
appear to be favored by observations. Several scenarios that 
enable the formation of low-mass stars in metal-free environ- 
ments have been proposed, among which the bi-modal IMF 
(Nakamura & Umemura 1999, 2002), enhanced HD-cooling 
in relic HII regions (Uchara & Irmtsuka 2000; Nakamura & 
Umemura 2002), the influence of a strong UV radiation field 
due to the vicinity of a very massive star (Omukai & Yoshii 
2003), fragmentation in shock-compressed shells induced by 



the first very massive SN explosions (Mackey, Bromm & 
Hcrntjuist 2003; Salvaterra, Ferrara & Schneider 2004). 

However, all these processes either require very special 
conditions, generally related to the existence of a previous 
generation of very massive stars, or have failed to be realized 
in realistic numerical simulations. In spite of this, observa- 
tions seem to indicate that either radiative feedback effects 
during the accretion phase on protostellar cores and/or a 
combination of the above mentioned processes can lead to 
the formation of Pop III stars with masses ^ IOOA/q. Ac- 
cording to the results of our analysis, the latter contribution 
to the Pop HI stellar IMF must be dominant at least at red- 
shifts z « 10. 

If this scenario is correct, relics of the low-mass Pop III 
star formation mode should be still present in our Galaxy, 
eventually contributing to the lowermost end of the metal- 
licity distribution function (Salvadori, Schneider & Ferrara 
in prep) . 
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